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Abstract 

Abnormal regulation of cell migration and altered rearrangement of cytoskeleton are characteristic of metastatic cells. The 
first described suppressor of metastatic processes is NI\/123-H1, which displays NDPK (nucleoside-diphosphate kinase) 
activity. To better understand the role of nm23 genes in cell migration, we investigated the function of NDK-1, the sole 
Caenorhabditis elegans homolog of group I NDPKs in distal tip cell (DTC) migration. Dorsal phase of DTC migration is 
regulated by integrin mediated signaling. We find that ndk-1 loss of function mutants show defects in this phase. Epistasis 
analysis using mutants of the ot-integrin ina-1 and the downstream functioning motility-promoting signaling module 
(referred to as CED-10 pathway) placed NDK-1 downstream of CED-lO/Rac. As DTC migration and engulfment of apoptotic 
corpses are analogous processes, both partially regulated by the CED-10 pathway, we investigated defects of apoptosis in 
ndk-1 mutants. Embryos and germ cells defective for NDK-1 showed an accumulation of apoptotic cell corpses. 
Furthermore, NDK-1 ::GFP is expressed in gonadal sheath cells, specialized cells for engulfment and clearence of apoptotic 
corpses in germ line, which indicates a role for NDK-1 in apoptotic corpse removal. In addition to the CED-10 pathway, 
engulfment in the worm is also mediated by the CED-1 pathway. abl-l/AbI and abl-l/Abi, which function in parallel to both 
CED-lO/CED-1 pathways, also regulate engulfment and DTC migration. ndk-H-);abi-l (-) double mutant embryos display an 
additive phenotype (e. g. enhanced number of apoptotic corpses) which suggests that ndk-1 acts in parallel to abi-1. Corpse 
number in ndk-1 {-);ced-10{-) double mutants, however, is similar to ced-10{-) single mutants, suggesting that ndk-1 acts 
downstream of ced-10 during engulfment. In addition, NDK-1 shows a genetic interaction with DYN-l/dynamin, a 
downstream component of the CED-1 pathway. In summary, we propose that NDK-l/NDPK might represent a converging 
point of CED-10 and CED-1 pathways in the process of cytoskeleton rearrangement. 
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introduction 

The human nm23 [nme] gene family consists of ten members 
named after the first identified metastasis suppressor nm23-Hl [non 
metastatic clone 23). The metastasis suppressor function has been 
extensively corroborated using metastatic cell lines (melanoma, 
breast-, colon-, hepato- and oral squamous cell carcinoma) where, 
for the most part, overexpression of NM23-H1 was associated with 
reduced cell motility (reviewed in [1]). Proteins encoded by the 
nm23 family are classified into two groups. Isoforms of group I 
(NM23-H1-NM23-H4) possess nucleoside diphosphate kinase 
activity and are highly conserved in eukaryotes from yeast to 
mammals [2] . Beyond their nucleoside diphosphate kinase activity, 
additional molecular functions are associated with NDPKs such as 



histidine-dependent protein kinase activity [3-4], 3 '-5' exo nucle- 
ase action [5-6] , DNase activity in caspase-independent apoptosis 
[7] and transcriptional regulation [8] . Together, group I members 
display essential functions; both up- and down-regulation can 
disrupt growth and/or differentiation [9; 10]. 

The most extensive set of studies analyzing group I members' 
role in cell motility and migration have utilized Dtmophita. awd, the 
fly orthologue of nm23-Hl / H2 is a negative regulator of migrating 
tracheal and border cells via modulating endocytosis of different 
receptors, such as platelet-derived growth factor receptor 
(PDGFR) /vascular endothelial growth factor receptor (VEGFR) 
[11] and fibroblast growth factor receptor (FGFR) [12]. In the 
process which affects the level of FGFRs Awd functions together 
with the dynamin/Shibire in endocytosis as a putative GTP 
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supplier for the GTPase [9]. Although no physical association of 
Awd and Shibire could be demonstrated in Drosophila, in rat a 
direct interaction was detected between NDPK and dynamin I by 

in vitro pulldown and coimmunoprecipitation [13]. Independent 
studies using Dictyostelium also confirm links to light-dependent, 
vectorial cell migration and cell nutrition through different forms 
of endocytosis [10]. 

Caenorhabditis ekgans serves as a particularly amenable model to 
investigate the process of cell migration. The nematodes are 
transparent and have simple anatomy making it possible to follow 
the migration of individual cells in the living animal throughout 
development. Well studied migrating cell types of C. ekgans include 
sex myoblasts (SM), two Q, neuroblasts (QL and QR) and their 
descendants, and distal tip cells (DTCs) or the gonadal leader cells 
[14-16]. In C. ekgans, we identified and described a single group I 
NDPK ortholog, NDK-1, which shows 85% and 86% similarity to 
NM23-H1 and H2, respectively [17]. In the current study we 
examine the role of NDK-1 in DTC migration and apoptosis. 

The tightly regulated migrator)- path of DTCs provides an 
elegant system to explore how cell migration can be guided within 
the spatial and temporal context of the organism. Distal tip cells 
are specialized leader cells and are responsible for gonad 
morphogenesis via their migration in C. ekgans [16]. During the 
four stages of larval development (LI, L2, L3, and L4), DTCs 
migrate in response to attractive and repulsive cues to properly 
form two mirror image U-shaped gonad arms (Figure lO). 
During development, somatic cells dying by apoptosis are engulfed 
by the neighboring cells as there are no speciedized engulfing cells. 
However, during germ cell death, which occurs as a part of the 
oogenesis program, at least half of all oogenic germ cells [17] are 
eliminated by apoptosis and are engulfed by gonadal sheath cells 
which surround the germ cells [18]. 

DTC migration and engulfment of apoptotic corpses are 
analogous processes in C. ekgans {Figure 1, [19-22]). Both require 
cytoskeletal rearrangements and membrane trafiGcking/recruit- 
ment [Figure lA, B) sharing genes regulating both processes 
{Figure ID, [22-23]). Two major partially redundant pathways 
regulate engulfment. The CED-lO/Rac pathway consists of unc- 
73/Trio, ced-2/CrkII, mig-2/RhoG, ced-5/DOCK180, ced-12/ELMO 
and ced-lO/Rac genes acting downstream of the alpha integrin 
receptor ina-1, while the CED-1 pathway includes ced-7/ABCl, the 
ced-1 /MEGFIO KCtpior:, ced-6/GULP, and dyn-1/ Dynamin [22-24]. 
The CED-1 0/Rac pathway controls both engulfment and the 
movement of DTCs by rearrangeing the cytoskeleton of the 
engulfing and migrating cells. The CED-1 pathway is involved 
only in engulfment where it recruits membranes to extend the 
surface of the engulfing cell. It has been suggested that a third, 
distinct pathway consisiting of abl-1 and abi-1 also influences both 
DTC migration and engulfment in parallel to the CED-10 Rac 
and CED-1 pathways [Figure ID) [22]. 

In this study we demonstrate that NDK-1 is required for normal 
DTC migration and engulfment of apoptotic corpses. We show 
that ndk-l influences both processes via common genes, acts 
downstream of ced-10 (ceU death abnormality) /-ftac and in parallel 
to abi-1 (Abl interactor)/AB/; and additionally shows a genetic 
interaction with dyn-1 /dynamin. Thus, NDK-1 affects the rear- 
rangement of cytoskeleton in both DTC migration and apoptotic 
engulfment. We also show that NDK-1 functions similar to its 
human counterparts in cell migration, as it inhibits the migratory 
potential of invasive breast adenocarcinoma cells. Our results 
might help to better understand the function of nm23 genes in 
metastasis. 



Results 

C elegans FLAG::NDK-1 reduced the motility of MDA-MB- 
2311 cells 

Our group is investigating the function of nucleoside diphos- 
phate kinases (NDPKs) in the model organism C. ekgans. NDK-1 is 
the single group I NDPK homolog of the worm and shows high 
sequence similarity to NM23-H1 and H2 [25]. It is known that 
nm23 genes regulate cell migration [26]. For example overexpres- 
sion of NM23-H1 and its sponge ortholog both reduced the 
migratory and invasive potential of CAL27 (oral squamous 
carcinoma of the tongue) cells [27]. Based on the high sequence 
similarity one might expect that the C. ekgans homolog of NM23- 
H1/H2 is also able to act likewise. Therefore we investigated the 
effect of NDK-1 exerted on the cell migration capacity of the 
breast adenocarcinoma MDA-MB-231T ceU line. MDA-MB- 
23 IT cells are far more migratory than CAL27 cells, and the 
influence of NM23-H1 is much more obvious in these cells. Stably 
transfected MDA-MB-231T cells overexpressing FLAG::NDK-1, 
FLAG::NM23-H1 and MYC-NM23-H2 {Figure 2A,B,C) were 
used for migration assay. Three independent experiments demon- 
strated that NDK-1 clones CEl and CE2 both diminished the 
migratory potential of MDA-MB-231T cells in a similar manner to 
overexpressing NM23-H1 (clones HAl and HA2) or NM23-H2 
(HBl and HB2) clones {Figure 2D). The suppression of migratory 
potential reached or exceeded 50% in almost all clones 
overexpressing either the worm or a human NM23 homolog 
compared to MDA-MB-231T control clones (clones Kl and K2). 

ndk-l(lf) mutants show incomplete migration of DTCs 

We recentiy characterized the worm ortholog of group I 
NDPKs, ndk-l, and identified the strong loss of function allele 
ok314. ok314 is a 1 157 bp-long deletion, which removes the entire 
ndk-l ORE, as well as upstream and downstream regulatory 

sequences [25]. 

ndk-l{ok3 1 4) mutants show a Pvl (protruding vuh'a) phenotype 
and they are sterile due to germ cell arrest in mitotic phase [25]. 
Morphological studies using Nomarski optics revealed that ndk- 
l{ok314) animals have insuflicientiy elongated gonad arms 
suggesting defects in DTC migration. Detailed analysis 
{Figure 3G,H) of 318 gonad arms showed that in tlu- majority 
of ok314 mutants (60.3%) DTCs turned to the proper direction 
and side however their migration is stopped prior to reaching the 
vulva. In 14.5% of the cases DTCs turned back to the vulva but in 
the ventral (instead of the dorsal) side or initiated ventral migration 
but subsequently vectored to the dorsal side before the turn (7.2%). 
Other defects manifested low penetrance (wrong direction, lack of 
the turn, wandering, bizzare twists) and we observed normal DTC 
migration in only 4.7% of the animals. Altogether we conclude 
that the prominent phenotype of ok314 mutants is incomplete 
migration of DTCs. 

In order to determine whether NDK-1 is expressed in DTCs, 
we studied NDK-1 expression in transgenic strains TTV2 and 
TTV3. Both were generated using the same translational construct 
but by distinct methods. In case of TTV2, the strain was created 
by ballistic bombardment, which resulted in extrachromosomal 
arrays of the transgene and therefore genetic mosaicism. For 
TTV3 the MosSCI method was appUed to integrate a single copy 
of the transgene into the genome. The latter method allowed the 
rescue of Ndk-l mutant phenotypes. Among TTV2 and TTV3 
animals we observed expression in DTCs from L3 {Figure 3 A, B) 
to L4 {Figure 3C,D) stage while in the isolated gonads of TTV3 
animals DTC expression was also detected in adults {Figure 3E,F). 
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Figure 1. Site of action of NDK-1 in the analogous processes of engulfment and DTC migration. A, B: Cell corpse engulfment and DTC 
migration are similar processes. In each case, the surface membrane of a cell (black) extends along the surface of another cell (hatched). The small 
arrows near the black cells indicate the directions of cell-surface extension. B: Only the relevant parts of body muscles are shown [1 9]. C: Schematic 
review of DTC migration (based on [21]). DTCs are located on the distal edges of the gonad primordium and start to migrate in L2. They migrate 
along the ventral surface (dashed line) of the hermaphrodite in L2 (first or ventral phase). Then they turn to the dorsal side during L3 (second or 
ventral to dorsal phase). A second turn redirects migration along the dorsal surface toward the center of the nematode during L4 (third or dorsal 
phase). The end of the migration is dorsal to the vulva, resulting in the mirror image U-shaped gonad of the adult. The developmental stage is 
indicated at right of each diagram. D: Signaling pathways in engulfment and DTC migration (based on [21-22]. Common genes are blue, green 
colour indicates the factors involved only in DTC migration, genes in purple boxes play a role only in engulfment. We suggest that NDK-1 /NIVI23 acts 
downstream of CED-10/Rac in the processes of DTC migration and engulfment of apoptotic corpses. NDK-1 shows a genetic interaction with DYN-1/ 
Dynamin. 

doi:1 0.1 371/journal.pone.0092687.g001 



As ndk-1 is expressed in tlie distal tip cells in difTerent 
developmental stages and knockout of ndk-1 results in difTerent 
DTC migration defects {Figure 3G,H), we hypothesized that ndk- 
1 functions in the process of DTC migration. 

ndk-1 {-) reduces the extra turn phenotype of CED-10 
pathway mutants while abi-1 restores the Ndk-1 DTC 
migration phenotype 

Integrins are heterodimeric receptors consisting one a and one 
P subunit. They provide connections between the extracellular 
matrix and the actin cytoskeleton regulating cellular shape, 
motility and cell cycle. Integrin mediated signaling is well known 



to be involved in the process of DTC migration [21]. Interestingly, 
both NM23-H1 and H2 have been linked to integrins [28-30]. 

The C. elegans genome encodes two a [ina-l: integrin alpha, pat- 
2: paralysed arrest two-fold embryo stage) and one fi (fiat-3) 
subunits [21]. AH these genes are expressed and function in DTCs 
[21;31-32]. Genes encoding C. elegans integrins are essential, null 
mutations in any of them cause embryonic lethality [21,31]. Since 
hypomorph integrin mutants are viable, we applied these mutant 
alleles for DTC analysis. Hypomorph ina-1 mutants show an 
abnormal migration (i. e. the migratory path of DTCs is often 
longer compared to wild-type and extra turns occur due to 
pathfmding defects, see also definition in methods) phenotype 
because DTCs do not stop prior to reproduction [2 1] . Silencing of 
pat-2 by RNAi resulted in dorsal pathfinding defects [21]. An 
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Figure 2. Western blot analysis and migration assay of transfected MDA-MB-231T ceils. IVIDA-MB-231T cells were stably transfected with, 
pcDNAS (K1 and K2), pcDNA3/FLAG-nm23-H; (HA1 and HA2), pcDNA3/FLAG-nd/<-/ (CE1 and CE2) and pcDNA3/MYC-nm23-H2 (HB1 and HB2). A: 
Western blot with anti-ot-tubulin antibodies (loading control). B: Western blot with anti-FLAG- antibodies, visible band in HA1, HA2, CE1 and CE2 
proves stable overexpression of introduced transgenes. C: Western blot with anti-MYC- antibodies, visible band in HBl and HB2 (overexpression of 
NIVI23-H2). D: Migration assay. MDA-MB-231T cells stably transfected with one of the following constructs: pcDNA3 (K1 and K2), pcDNA3FLAG/nm23- 
H1 (HAl and HA2), pcDNA3FLAG/nd/(-? (CE1 and CE2) and pcDNA3/MYC-nm23-H2 (HBl and HB2) were tested for migration potential. The cells were 
stained with crystal violet and counted (the number of migrated cells were counted in four representative microscopic fields per each clone). The CE1 
and CE2 clones as well as HAl, HA2, HBl and HB2 exhibited significantly diminished migration potential compared to control (K1 and K2) clones 
(Student's t-test, p<0.05). The results are presented as an absolute number of migrated cells in 4 representative fields for every clone (±SD). 
doi:1 0.1 371 /journal.pone.0092687.g002 



incomplete migration with an enlarged and blunt end of the gonad 
arm can be observed with DTC-specific expression of a dominant- 
negative pat-3 transgene or by pat-3{RNAi) [21,31-32]. 

vab-3 (variable abnormal morphology) encodes a homeodomain 
protein, the C. elegans orthologue of Pax6. vab'3/Pax6 transcrip- 
tionally regulates both ot-integrins in different manners: it 
downregulates ina-1 [Figure ID) to cause the cessation of DTC 
migration and activates pat-2 expression at L3 stage which is 
necessary for normal dorsal pathfinding [21]. vah-3 reduction of 
function mutants show an overmigration phenotype similar to ina- 
lirfj mutants. 



The Rac GTPases CED-10 and MIG-2 (abnormal cell 
mig ration) have been shown to act downstream of INA- 1 during 
the migration of DTC s [21]. Additionally, ced-10 was found to 
function together with ced-2 and ced-5 in DTC migration (reviewed 
in [16]). Moreover, in the last decade, an entire pathway (often 
called as CED-10 pathway) was built using also the above 
mentioned genes [Figure ID). [19,21-22]. Mutations in CED-10 
pathway components manifest characteristically extra turns due to 
pathfinding defects. Furthermore, these genes - similar to integrin 
signaling components - function during the dorsal (or third) phase 
(Figure IC) of DTC migration [16]. 
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Figure 3. NDK-1 functions in DTCs. NDK-1 ::GFP shows expression in distal tip cells (DTCs) in TVJ2 (B, D) and TTV3 (F) translational reporter lines 
in L3 (A, B), L4 (C, D) larvae and adults (E, F and I). A, C, and E are the corresponding DIC images of B, D, and F, respectively. DTC locations are 
indicated with arrows. We note that on panel D GFP expression in DTC (marked by an arrow) is overshadowed by intense expression of 
autofluorescent granules in the intestine. E, F and I: gonads were isolated from adult animals. H: DIC image of an adult wild-type (N2) gonad arm. G: 
Variation and distribution of DTC migration phenotypes observed in ndk-l(ok314) mutants (see details in the main text). The migratory path of DTCs 
is marked by dashed lines. I: In TTV3 lines carrying the integrated transgene, NDK-1 ::GFP expression is also observed in gonadal sheath cells. DIC, 
fluorescent and merged images derive from the gonadal loop region of a transgenic animal. On the DIC panel white arrows show dying cells, which 
are surrounded by sheath cells strongly expressing NDK-1 ::GFP (white arrows on the fluorescent and merged panels). 
doi:1 0.1 371/journal.pone.0092687.g003 



ABL-1 and ABI-1 function in a distinct, recently discovered 
cascade, acting parallel to the CED-10 pathway in DTC migration 
{Figure ID, [22]). abl-l{okl71) and abl-l{nl963) alleles manifest no 
obvious gonadal phenotype [22] but the DTC migration defects of 
ced-2, ced-5, ced~12 and ced-10 mutants were suppressed in abl-1 
mutant backgrounds. ABL-1 inhibits ABI-1 in C. elegam, abi- 
l[tm494) mutation itself or abi-l[KNA{) caused only weak defects in 
DTC migration [22], however abi-l{RNAi) enhanced the extra 
turns oi ced-5[nl812) and ced-12{n3261) [22]. We analyzed DTC 



migration in abi-l[ok640) mutants and observed 26% of abnormal 
gonad arms in which the overshoot phenotype was the most 
frequent (13%) (Figure 4 A, B). 

NDK-1 expression was detected in DTCs from L3 larval stage 
onwards, and lack of NDK- 1 activity resulted in early cessation of 
migration in the dorsal phase in most of the cases [Figure 3G). 
Based on these facts we investigated the connection between ndk-1 
and integrin signaling genes during the process of DTC migration. 
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Figure 4. nt/^- 7 functions lilieiy downstream of ced-10'm tKie process of DTC migration and shows a genetic interaction with abi-1. A: 

Epistasis analysis using mutants of the homeodomain transcription factor vab-3, the a-integrin ina-1 and the downstream functioning CED-10 
pathway genes suggests that NDK-1 acts downstream of CED-10/Rac. Loss of ndk-1 reduced the extra turn phenotype of vab-3, ina-1, unc-73, mig-2, 
ced-5, ced-12 and ced-10 mutants (B) and mostly resulted in Ndk-1-like gonad (B). Examining double mutants of ndk-1 and the parallel pathway 
functioning in DTC migration, we found that abl-1(ok171) did not influence the reduced migration of ndk-1 mutants. However, the abi-1(ok640) 
mutation (B) was able to partially restore the DTC migration defects of ndk-1{ok314] animals (B). The different DTC migration categories are described 
in iVlaterials and Methods. Statistical analysis used Fisher's exact test. The difference in the distribution of migration phenotypes was significant 
(Fisher's exact test, psO.05) in all pairs of single and double mutants (with the exception of ndk-1(RNAi) and vab-3(e1796);ndk-1{RNAi) mutants 
(p>0.05)). 

doi:1 0.1 371 /journal.pone.0092687.g004 



We performed an epistasis analysis (a genetic tool suitable for 
ordering genes into pathways) using the ok314 mutation or RNAi 
mediated silencing of ndk-1 and the following loss of function or 
reduction of function mutants: vab-3(el796), ina-1 (gm39), ced- 
10(nl993), mig-2(lql3), unc-73(e936), ced-5(nl812), ced-12(n3261), 
abl-l(okl71)/ASLl, abi-l[ok640). Except ahl-l[okl71) and ahi- 
l{ok640), all these mutants display large gonads with extra turns 
(Figure 4A). 

To examine the effect of the absence of NDK- 1 in the above 
mentioned mutant backgrounds, we either depleted the function of 



NDK-1 by RNAi or generated double mutants. In the cases oivah- 
3 (el 796), ina-1 (gm39), unc-73 (e936), ced-5(nl812) and ced-12{n3261) 
RNAi was applied to knockdown NDK-1 function in the 
corresponding mutant background. In the cases of mig-2(lql3), 
ced-10(nl993), ahl-1 (okl 71) and ahi-l(ok640) double mutants were 
created by introducing the ok314 mutant allele of ndk-1 into the 
appropriate single mutants. 

Loss of ndk-1 suppressed the extra turns of mb-3, ina-1, unc-73, 
mig-2, ced-5, ced-12 and ced-10 mutants and mostiy resulted in Ndk- 
1 -like gonad (reduced gonad with fewer extra turns) [Figure 4 A, 
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B). Although, as mentioned above, abl-l[okl71) was able to 
suppress the phenotype of the c«rf mutants [22], it did not influence 
the reduced migration of ndk-l mutants. Unlike abl-l{okl71), ahi~ 
Uok640) mutation was able to restore the DTC migration defects 
of ndk-l(ok314) animals in of the cases {Figure 4 A, B), 

suggesting that ndk-l might either act upstream of abi-1 or loss of 
ABI-1 compensates for the absence of NDK-1, indicating parallel 
functions of the two genes. 

Based on these results we suggest that ndk-l might act 
downstream of or in parallel to ced-lO/Rac and upstream of or 
in parallel to abi-1 /Ahi [Figure ID) in the process of DTC 
migration. Our data raise the possibility that ndk-l acts 
downstream of ced-10 and in parallel to abi-1. Although this 
tendency could be hypothesized on the basis of the distribution of 
migration phenotypes, it could not be confirmed by statistical 
analyses {Figure 4A). Thus, the epistatic relationship of NDK-1 
and the two parallel pathways cannot be deduced unambiguously, 
as DTC migration phenotypes are difflcult to quantify. To specify 
NDK-l's site of action more precisely, we decided to analyze the 
role of NDK-1 in apoptosis, where phenotypes, (e.g. number of 
apoptotic corpses, so called Ced phenotype) can be quantified with 
precision. 

ndk-l {-);abi-l {-) double mutants show an additive Ced 
phenotype, suggesting that NDK-1 acts downstream of 

CED-10, in parallel to ABI-1 

Since CED-lO/Rac signaling controls both DTC migration and 
engulfment, and our data derived from the DTC analysis 
suggested that NDK-1 acts downstream of or in parallel to 
CED-10, we hypothesized that NDK-1 also functions in the 
engulfment phase of apoptotis. In somatic tissues of the worm, 
elimination of apoptotic corpses occurs by the non-specialized 
neighbouring cells, however cell corpses in the germline are all 
engulfed by gonadal sheath cells [17-18]. Transgenic worms of the 
strain TTV3 show that NDK-1 ::GFP is expressed in these somatic 
gonadal cells around dying germ cells (Figure 31), which further 
supports the potential role of NDK-1 in engulfment. To address 
this issue, first we examined the number of germ cell corpses in 
ndk-l mutants. Among germ cells, only the developing oocytes die 
via apoptosis when they exit the pachytene phase at the gonadal 
loop region [17-18]. 

We used RNAi because germ cells of ndk-l{ok314) loss of 
function mutants arrest at the mitotic phase before reaching 
pachytene (our unpublished data), thus engulfment defects cannot 
be investigated in this mutant background. We analyzed the 
gonads of ndk-l{RNAij worms using DIG optics and also marked 
the early apoptotic corpses with CED-1::GFP [24]. \Sk- obsi^rved 
an excessive level of GFP-labeled apoptotic corpses in ndk-l{RNAi) 
treated worms (6.3; N = 87) compared to the control animals (2.9; 
N = 65) {Figure 5A-E). 

During C. elegans development, 1 3 1 somatic cells are eliminated 
by apoptosis [33] and most of them (109) die during mid- 
embyrogenesis [34] . Expecting that the lack of NDK- 1 results in 
an accumulation of cell corpses during embryogenesis as well, we 
counted apoptotic corpses by DIG optics in ndk-l(ok314) and in 
wild type embryos around the comma stage. We observed a near 
doubling of corpses in ndk-l mutant embr\'os (20.1) compared to 
wild type (10.9), indicating that NDK-1 plays a role in apoptosis 
{Figure 5F,G,L and panels A,D in Figure SI). 

Next, we examined NDK- 1 's function in the engulfment process 
during embryogenesis. We analyzed apoptotic corpses in comma 
stage embryos and in embryos slightly after comma stage (close to 
1.5 fold stage) harboring mutations in ndk-l and the key positioned 
genes abi-1 and ced-10 (downstream components of the CED- 10/ 



Rac and ABL-l/ABI-1 parallel pathways). First we counted 
apoptotic corpses in the homozygous mutant progeny of ndk- 
I(ok314)/+;ced-10(nl993) animals and found that double mutant 
comma stage embryos contained on average 25.2 apoptotic 
corpses {Figure 5 J, G, K, L, M and panels A,C,F in Figure 
SI) which is comparable with ced-10{nl993) single mutants; 24.8. 
The majority of the 25 examined ndk-l{ok314)\ced-10{nl993) 
double mutant embryos displayed cell corpses between 22-27, 
reminiscent of ced-10{-) single mutants, however around 10% of the 
doubles showed an increased number of cell corpses (e. g. in two 
embryos 38 corpses were counted) {Figure 5 L, M). We note that 
for analysing the genetic interaction between ced-10 and ndk-l, the 
ced-10 allele nl993 was used, which is considered as a partial loss- 
of-function allele, not a null, because ced-10 nuU mutations are 
maternal effect lethal [35]. 

Next, we examined the genetic interaction of NDK-1 and ABI-1 
in apoptosis. abi-l{tm494) or abi-l{RMAi) alone caused only weak 
defects in engulfment [22]. Consistently, we detected slightly more 
apoptotic corpses in abi-l{ok640) embryos (18.7) compared to wild 
type (10.9) {Figure 5 F, H, L, M). 

Analyzing the homozygous progeny of ndk-l (ok3 14)/ +;abi- 
l(ok(>40) heterozygotes revealed that ndk-l{ok3 1 4);abi-l{ok640) 
double mutants displayed 36.7 cell corpses on average {Figure J 
F, G, H, I, L, M, panels B,D,E in Figure ST). Thus, 
introducing abi-l{ok640) mutation into ndk-l{ok314) background we 
noticed an additive phenotype, an increased level of apoptotic 
corpses {Figure 5 1, L, M). This additive phenotype suggests that 
ndk-l acts in parallel to abi-1. Namely, if two engulfment genes act 
in the same linear pathway, the phenotype of the double mutants 
(number of corpses) should not be more severe than that of the 
stronger single (nuU) mutant. If the two genes function in parallel 
pathways a (nuU) mutation in one of them should enhance 
significantiy the phenots-pe caused by the other, resulted in an 
additive phenotype [24,36] . Therefore our results suggest that in 
the process of engulfment ndk-l acts in parallel to abi-1 and 
downstream oi ced-10 {Figure ID). 

nclk-Hok314) embryos show a phenotype characteristic 
for dyn-li-) mutants: late embryonic lethality with 
persistent cell corpses 

NDPK/AWD is known as a potential GTP suppher of the large 
GTPase, Dynamin/Shibire in several model systems [11-13,37]. 
The C. elegans ortholog ot dynamin/shibire, dyn-1 acts in the GED-1 
engulfment pathway downstream of ced-6 ([36]; Figure ID). 
DYN-1 is proposed to organize vesicle transport: 1. to the 
phagocytic cups for extending pseudopods during engulfment and 
2. to phagosomes for apoptotic cell degradation [36]. 

dyn-1 loss of function mutants are embryonic lethal [36,38] due 
to endocytosis defects. Beyond this main phenotype, defective 
DYN-1 causes failure in engulfment [36]. Many djn-1 mutant 
alleles bearing missense mutations in the GTPase domain result in 
embryonic lethality at later embryonic stages and a strong Ced 
phenotype. These embryos do not move inside the eggshell in 
contrast to the other same-aged embryos. The combination of 
persistent cell corpses and late-embryonic lethality was defined as a 
new phenotypical class [24,36] . 

Previously wi; characterized Ndk-l mutant phenotypes and 
showed that in the progeny of ndk-l (ok3 1 4) /+ heterozygotes 12% 
of ndk-l {ok3 14) homozygotes die as embryos [25]. Further analysis 
revealed that these dying ndk-l{ok314) embryos also show the Dyn- 
1-like late embryonic lethality phenotype with persistent cell 
corpses {Figure 6B). 

Next, to see whether ndk- 1 and dyn-1 also interact in the worm, 
we decided to generate ndk-l{-);d)>n-l{-) double mutants. For 
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M 



genotype 


n corpses ± sd 


n 


wild type 


1^10,9 ±0,7 


10 


ndk-l(ok314) 


20,1 ± 1,7 


15 


abi-l(ok640) ^^^^^ 


^ 18,7 ± 2,0 


15 


ndk-l(ok314);abi-l(ok640) 


36,7 ±3,2 


10 


ced-10(nl993) ^^^^^ 


B 24,8 ± 2,2 


20 


ndk-l(ok314);ced-10(nl993) 


25,2 ±4,0 


25 




Figure 5. ndk-l(ok314);abi-l(ok640\ double mutants show additive Ced phenotye. Pn^.jced-lr.gfp transgenic worms treated by ndk-HRNAi) 
(D) show an excess of apoptotic corpses in the germline compared to worms carrying the same transgene treated by control RNAi (B). A, C are 
corresponding DIC images of B, D respectively. Arrows indicate apoptotic germ cells. E: Increase of apoptotic germ cell death in ndk-1{RNAi) animals 
compared to the control, where pn^.jCed-V.-.gfp transgenic worms were grown on control RNAi (e.g. £ coti HT1 15(DE3) carrying an empty vector). F- 
K; Monitoring apoptotic corpses in wild-type embryos (F), ndk-H-) (G), abi-l{-) (H), ced-10(-) (J) single mutants and ndk-1 (-);abi-H-) (I), ndk-1 (-);ced-10(- 
) (K) double mutants using DIC optics. Embryos slightly before or at the comma stage were scored. Each panel shows two focal planes (F-K). 
Arrowheads indicate apoptotic corpses. Panel L shows a summary of apoptotic corpses scored in ndk-l{ok314), abi-l(ok640), ced-10(nl993) single 
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mutant and ndk-1 (okSl 4);abi-1 (ok640) and ndk-1 (ok31 4):ced-l 0(n 1 993) double mutant embryos. ndk-l(-);abi-l{-) double mutants display an enhanced 
Ced phenotype compared to single mutants, however ndk-1 (-);ced-10(-) doubles are reminiscent of ced-10(-] single mutants. Panel Wl is the graphic 
representation of panel L. p-values refer to comparisons of apoptotic cell corpse numbers between single and double engulfment mutants. * * * 
means p<0.001; n.s. means not significant (p = 0.627). 
doi:10.1371/journal.pone.0092687.g005 



crossing we used ndk-l{ok314) mutants and the viable kfSl allele of 
dyn-1, wherein we took the advantage of thermosensitivity. 

We performed this experiment at the restrictive temperature for 
kySl (25°C) {Figure 6A). 15.8% of ndk-1 {ok3 14) homozygotes died 
as embryos and ndk-l/-\- heterozygotes laid on the average 122.7 
eggs. We measured approximately 50% lethality of djn-l{ky5T) 
single mutants and 51% lethality in the Fl progeny of djn- 
l(kj51);ndk-l(ok314)/-\- animals, but the broodsize was decreased 
(43.9) in the latter category compared to dyn-1 single mutants (63). 
At the restrictive temperature we did not notice any Pvl, Ste 
(double mutants) in the Fl progeny. Putting the Fl animals to the 
permissive temperature (15°C) the Pvl, Ste phenotype appears in 
F2 generation. Therefore these observations suggest that the 
double mutant is not viable indicating that ndk-1 and dyn-1 interact 
genetically in the worm as well. 

Discussion 

Group I NDPKs are negative regulators of cell migration and 
motility. To investigate the function of the sole worm group I 
homolog, NDK- 1 in cell migration, we overexpressed NDK- 1 and 
human NM23-H1 and H2 in MDA-MB231T, an invasive breast 
carcinoma cell line. We found that NDK-1, similar to its human 
counterparts, dramatically suppressed the migratory potential of 
MDA-MB231T ceUs. This result shows that NDK-1 bears an 
evolutionary conserved function in cell migration and that C. 

A 



ekgans serves as a tractable model system to monitor the effects of 
the NDPK gene family. 

Next, we studied the role of ndk-1 in the process of the migration 
of distal tip cells (DTCs) which are responsible for leading gonad 
morphogenesis. During larval development distal tip cells guide 
the migrating gonad arms which finally form two U-shaped tubes 
[16]. Studying different transgenic strains, we detected NDK- 
1::GFP expression in DTCs. Detailed analysis of ndk-1 loss of 
function mutants showed different DTC migration defects, but 
predominantly, incompletely elongated gonad arms were the 
consistent feature. The majority oi ok314 mutants showed a third 
(dorsal) phase migratory defect, as a consequence of a premature 
stop of DTCs, after the turn, that was not present in wild-type 
worms. 

Thus, overexpression of NDK-1 inhibited cell migration and 
reduced motility of metastatic breast adenocarcinoma cells, but 
loss of NDK-1 function in the worm predominantly caused 
incomplete gonadal migration and not excess migration. These 
data suggest that NDK-1 might have opposing functions 
depending on the cellular environment, e.g. NDK-1 might 
promote cell migration in the worm, but inhibit it in human cell 
lines. In support of this idea, recent work showed that NDPK 
affects Dictyostelium growth in opposing ways in difierent cellular 
environment [10]. To further investigate this issue, it wiU be 
important to overexpress NDK-1 in the worm and analyze the 
phenotypes in response to overexpression. 



Genotype of the percentage of embryonic percentage of Pvl, expected ratio of average brood 

hermaphrodite (P) lethality In Fl Ste progeny (Fl) Pvl, Ste in Fl (%) size of Fl 

ndk-l(ok314)/+ 15,8 9,2 25 122,7 

dyn-l(ky51) 50 ND ND 63 

ndk-l(ok314)/+;dyn-l(kySl) 51 0 9^ 43^9 




Figure 6. ndk-1(ok314l;dyn-1(ky51) double mutants are lethal. A: In the Fl progeny of nc//<- !(-)/+ heterozygotes only 9.2% homozygous Pvl, Ste 
adults can be observed instead of the expected 25%, since 15.8% of ndk-1 (ok314) homozygotes die as embryos at 25°C. At the restrictive temperature 
(25 C) 50% of dyn-l(ky51) single mutants die as embryos. At 25 C, in the Fl progeny of dyn-l(ky51);ndk-1(ok314)/-\- animals we got decreased brood 
size and we did not notice any Pvl, Ste animals, suggesting that the double mutants are not viable. B: 3-fold stage homozygous ndk-l(ok314) embryo 
shows late embryonic lethality with persistent cell corpses. Arrows indicate apoptotic corpses. 
doi:1 0.1 371/journal.pone.0092687.g006 
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Previously nm23 genes were linked to intcgrin signaling in 
human studies [28-30]. In C. elegans, integrins act in the dorsal/ 
third phase of DTC migration [16,21], where the majority of 
defects conferred by ndk-l(ok314) mutants are also observed. 
Consistently, NDK-1::GFP expression was detected in the third 
phase of migration (from L3 to adulthood). Based on these findings 
we performed an epistasis analysis with genes involved in CED- 
10/Rac signaling {pab-3/Pax6, ina-1 /alpha-integrin, unc-73 /Trio, mig- 
2/RhoG, ced-5/Dockl80, ced-12/Elmo, ced-lO/Rac) and with ahl-1/ 
ABL and ahi-l/ABI which were shown to act in parallel to the 
CED-IO pathway. All mutants of CED-IO pathway genes have 
large gonads with extra turns. Putative null mutant alleles of abl-1 
(okl71 and nl963) display no obvious gonad phenotype [22]. We 
observed 30% of abnormal gonad arms in abi-l{ok640] mutants in 
which the overshoot phenotype was the most frequent (14%). This 
result is consistent with previous data, e.g. ahi-1 mutations or abi- 
l(RNAi) results in overextensioned DTC migration [22]. Loss of 
ndk-l reduced the extra turn phenotype oivab-3, ina-1, unr-73, mig- 
2, ced-5, ced-12 and ced-10 mutants and mosdy resulted in Ndk-1- 
Uke gonad (similar to that of ndk-l[-) mutants). These results lead to 
the hypothesis that NDK-1 might act downstream of or in parallel 
to ced-10. Although ahl-l(()kl7 1) was able to suppress the extra turn 
phenotype of the mutants [22], it did not influentx tlu; rc-duced 
migration of ndk-l mutants. Unlike abl-l{okl71), abi-l(pk640) 
mutation was able to restore partially the DTC migration defects 
oindk-l{ok314) animals, showing a genetic interaction between ndk- 
1 and abi-1. 

Both distal tip cell migration and engulfment of apoptotic cell 
corpses require precise regulation of cytoskeleton to extend cell 
surfaces [19-20,22]. Although nm23 has been implicated in 
apoptosis [7,39-43] in these studies NM23's function was detected 
in the dying, not in the engulfing cells. In C. elegans, during somatic 
cell death neighboring cells engulf the nascent apoptotic corpses; 
however germ cell corpses are eliminated by specialized engulfing 
cells, the sheath cells [17]. We noticed increased numbers of 
apoptotic corpses in the germline of ndk-l{RNAi) animals and in 
ndk-l{ok314) embryos, furthermore, NDK-1 is expressed in sheath 
cells (this study) as well as in embryos [25]. These results suggested 
that NDK-1 plays a role in engulfment of apoptotic cells. We note 
that subsequent experiments are necessary to further confirm the 
role of NDK-1 in engulfment, such as investigating the persistance 
of apoptotic corpses in ndk-l{-) single mutants. 

Starting from the observation that ndk-l genetically interacts 
with abi-1 in DTC migration and functions likely downstream of 
ced-10, we hypothesized that ndk-l regulates engulfment in the 
same manner. We analyzed ndk-l (-), abi-1 (-) and ced-lO(-) single 
mutant embryos and ndk-l (-);ahi-l(-) and ndk-l (-);ced-10(-) double 
mutants. ced-10{-), ndk-l{-) and abi-l{-) mutants contain more 
apoptotic corpses than wild type animals. Moreover, we observed 
an increased number of corpses, an additive phenotype in ndk-l{- 
);abi-l[-) double mutant embryos, whereas ndk- 1 {-);ced- 1 0{-) double 
mutants displayed apoptotic corpses reminiscent to ced-10[-) single 
mutants. These results suggest that ndk-l acts in parallel to abi-1 
and downstream of ced-10 during apoptotic enguUinent. 

We also analyzed the putative interaction of ndk-l and dyn-1/ 
dynamin, which is a downstream factor of the CED-1 pathway 
involved in engulfment. The CED-1 pathway is thought to recruit 
membranes to extend the surface of the engulfing cell, and DYN- 1 
is proposed to organize \esicle transport during this process. 
Interaction between NDPKs and Dynamins in endocytosis is well 
known in other model systems [9,12-13,36,44]. ndk-l{-) and dyn-l{- 
) single mutants display a phenotypic similarity in the worm, the 
late-embryonic lethality with persistent cell corpses. Double 
mutant analysis revealed that ndk-l(ok314);dyn-l(k)/51) animals 



are not viable indicating a genetic interaction between the two 
genes in C. elegans as well. 

Togeth(;r, we propose that NDK-1 acts downstream of CED-10 
and in parallel to ABI-1 in the process of engulfment and DTC 
migration. In addition, ndk-l shows a genetic interaction with dyn- 
1/ dynamin, which raises the possibility that NDK-1 might represent 
a converging point of the CED-1 and CED-10 pathways. 
Interestingly, in 10% of ndk-l {-),ced-10{-) double mutant embryos 
we observed an increased number of cell corpses (note that in 90% 
of the double mutants the corps number is similar to ced-10{-) single 
mutants). The c()n\'(;rgc'ncc of tiu; two pathways would explain this 
partial additivity observed in ndk-l{-);ced-10{-) double mutants. 
There are also data showing that silencing of rfvw-i results in a 
DTC migration phenotype [38], although the CED-1 pathway is 
not considered to influence DTC migration. Earlier it was 
suggested that the CED-10 and CED-1 engulfment pathways 
converge at ced-10 [45], but later studies dealing with DYN-l's 
function presume that the two engulfment pathways do not 
converge [36,46]. Recendy, it was found that the CED-1 pathway 
not only regulates membrane recruitment during engulfment, but 
is also connected through CHC-l/clathrin heavy chain to F-actin 
assembly [46]. Thus, the CED-1 pathway - besides the CED-10 
pathway — is also linked to cytoskeleton remodeling. Earlier studies 
confirmed NDPKs' role in cell migration and endocytosis [12,26], 
processes, which are both related also to CED-1 and CED-10 
pathways. 

In summar)', our geni;tic analysis shows that ndk-l acts 
downstream of ced-10 in the regulation of engulfment and DTC 
migration and genetically interacts with dyn- 1 / dynamin. These data 
might contribute to our understanding of how NDK-1 exerts its 
function in the rearrangement of the cytoskeleton. Although 
NDPKs were already related to phagocytosis in different systems 
[10,47], this is the first time that we link NDK-l/NDPK to 
apoptotic engulfment. We show that impaired NDK-1 function 
cause's inefficient engulfment. In mammals it is known that failure 
of engulfment leads to inflammation [48] , moreover inflammation 
favors tumour progression in certain circumstances [49] . In further 
studies it would be interesting to overexpress or silence NDK-1 in 
specialized engulfing cells, such as the macrophage-monocyte 
system, to further investigate NDPKs' function in apoptotic 
engulfment. 

Materials and Methods 

Nematode strains and alleles 

C. elegans strains were cultivated at 20°C on NGM plates seeded 
with OP50 bacteria [50]. The N2 Bristol strain was used as the 
wild-type strain. The following transgenic or mutant strains were 
used: 

TTVl ndk-l(ok314) I/hT2[bU-4(e937) let-?(q782) qIs48](I;III); 
TTV2 eluExl[NDK-l::gfp; unc-1 19(+)];unc-l 19(ed3)III; TTV3 
eluSil[NDK-l::GFP+cb-unc-119(+)]III; LGI: MT11068 ced- 
12(n3261), CB3203 ced-l(el735), CB936 unc-73(e936), LGIII: 
RB829 B0336.6(ok640) corresponds to abi-1 (ok640), NG39 ina- 
l(gm39), LGIV: MT4434 ced-5(nl812), MT1522 ced-3(n717), 
MT5013 ced-10(nl993), LGV: MD701 bcIs39[P(lim-7)ced- 
l::GFP+lin-15(+)], LGX: CX51 dyn-l(ky51), CB3304 vab- 
3(el796), LE815 mig-2(lgl3), XRl abl-l(okl71). 

Information about alleles can be found at www.wormbase.org. 

RNA interference 

For RNAi through ingested dsRNAs, an Xhol digested 
fragment of the full-length cDNA corresponding to F25H2.5 
(ykll05e04) was cloned into pPD129.36, and the obtained 
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construct was transformed into HT115 (DE3) bacteria. RNAi 
experiments were performed at 25 °C essentially as described [51]. 

Quantitation of DTC defects. Adult animals were anaes- 
thetized and examined visually. Only completely visible gonads 
were analyzed. Gonadal length was defined as reduced when the 
gonad tip was before the vulva and as overshoot when the gonad 
was elongated past vulva. Extra turn category means that we 
observed more than one complete turn. DTC migration was 
scored as 'other' when the gonad was morphologically abnormal 
(wrong direction, lack of the turn, wondering, extra arm, bizarre 
twists). As animals affected by ndk-URNAi), sterile adults possessing 
a protruding vulva were picked. For statistical analysis DTC 
migration phenotypes were grouped into three classes: reduced 
migration, migration positioned to vulva and excess migration (in 
the last class overshoot and extra turn categories were merged, see 
Figure 3 A). Animals Ix^longing to the 'other' rategor)- were not 
included in the statistical analysis, as they display a broad 
phenot)'pic variance, cannot be grouped in any category, and 
their ratio does not reach 10% in any mutant population used, p- 
values of pairwise comparisons of DTC defect rates were 
calculated by Fisher's exact probability test. 

Quantitation of apoptotic cell corpses. Apoptotic corpses 
(as refractile discs) were scored in comma stage embryos directly 
by microscopy [33], using Nomarski/DIC (differential interference 
contrast) optics on Olympus BX51. /i-values for pairwise 
comparisons of apoptotic cell corpse numbers were calculated 
using the Student's t-test. 

Cloning of the FLAG::NDK-1 (pcDNA3FLAG-nm23-NDK) 
construct 

The insert of the construct was generated from the cDNA 
ykl 105e04 by PGR using the following primers: 5'-cag aag ate tat 
gga eta caa gga cga cga cga taa gat gag caa cac tga gag aac c-3' 
(which contains the flag sequence) and 5'-ata gtt tag egg ccg ctt tat 
teg tag ace cat gag ttg-3'. The insert was ligated into pCDNAS.l 
using BamHI and NotI sites. 

Preparation of stably transfected MDA-MB-231T clones 

MDA-MB-231T cells were transfected as follows: 4x10*' cells in 
DMEM supplemented with 10% FBS were subjected to electro- 
poration (250 V, 950 [xF) with 25 |xg of each of the following 
plasmids: pcDNASFLAG-wm^^-NDK, pcDNA3FLAG-«m25-//i, 
pcDNA3MYC-«m25-//2 and pcDNA3 as control. Post transfec- 
tion the cells were seeded on 4x100 mm plates and incubated in 
DMEM supplemented with 1 mg/ml geneticin (Sigma) until 
development of resistant colonies. Positive clones were screened by 
Western blotting, propagated and stored in liquid nitrogen until 
further usage. 

Western blotting 

For detection of stably transfected MDA-MB-231T clones 10*^ 
cells of each clone was collected, washed in PBS and sonicated. 

The samples (10 ng) were loaded on SDS-PAGE and electro- 
transfered to an Imobilon-PSQ_ membrane (Milipore). For 
detection of FLAG/NDK and FLAG/HI the anti-FLAG M2 
antibody (Sigma) was used while MYC/H2 was detected using 
anti-MYC antibody (Santa Cruz Biotechnology). Anti-a-tubuhn 
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antibody (Calbiochem) was used as loading control. After 
application of appropriate secondary antibodies the protein bands 
were visualized using Western Lightning Plus-ECL (PerkinElmer, 
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Supporting Information 

Figure SI Monitoring apoptotic corpses in embryos 
slightly after the comma stage in different mutant 

backgrounds. A-F: Monitoring apoptotic corpses in wild-type 
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Each panel shows two focal planes (A-F). Arrowheads indicate 
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